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Even with Incrcasocl recognition since the 1960 's of wetland jjfnpor- 
tancc, grov/ing world population and food demand have resulted in exten- 
sive v/ctland drainage. Since World VJar II nearly 25 ^^ of the v/etlands in 
the Prairie Pothole Region of the Dakotas cind Minnesota has been drained. 
Sound federal, state, and local management decisions on acquisition, pres 
ervation or drainage of v/etlands must be based on sound inventory data. 
Currently, the U.S. Fish and V/ildlife Service is involved in a National 
Wetland Inventory to satisfy these data needs at the federal level. Like 
wise, various state emd local agencies are undertaking, or have completed 
v;eltand inventories for a range of purposes. For example, in Minnesota 
v/etland inventories serve as baseline data to evaluate site modification 
permits, non-point pollution sources, and generally to execute the public 
water la\/. 

In short, accurate classification and mapping of v;etland cover types 
is considered essential to the implementation of a nimiber of federal, 
state, and local legislative mandates and resource management programs. 

At the some time, the task of max^ping v;etlands on a statewide basis is 
challenging in terms of cost, time and manpower. On a national level, 
the task is monumental. 

The general objective of the federal v/etland inventory is to clas- 
sify and locate the nation's wetlands as accurately as possible on 
1:100,000 (or in some cases 1:24,000) scale mai^s (Cov/ardin £t , 1977). 
To accomplish this, two xx)ssiblG reraotc sensing data bases suggest them- 
selves: (a) high altitude, small-scale aerial p>hotograx')hy and/or (b) 

LAMDSAT. Of the tv/o , cierial photogrciphy apx^ears to be the favored choice 
from the standjpoint of spatial resolution - but x^hotography is not o„lv;ays 
available ajid, when available, it not always adequate for the purpose. 


Ncithor tdolmicxuc, LAUDSAa’ in particular, had boon oufficiontly toafberl in 
arcar, reprosontativo of Minnesota i>xiox to thiu study. Thus, this study 
was proposed to test the applicability of automated processing of liANDSAT 
multispectral scanner (MSS) in the Minnesota regional contoxt. 


STUDY OBJECTIVES 


Initially, the sole objective of this study was to tost, under local 
conditions, the capability of Landsat data aricilysis techniques '• o xx>si- 
tion adequately and classify v^otlands in support of the Ncitior.au Wetlands 
Inventory. This involved a comparison of single-date vs. double-date 
sets as well as a comparison of data sots preprocessed with a coarse vs. 
a precision geometric correction. The test area chosen for these compar- 
isons consisted of five contiguous 7.5-ininute quadrangles located v/itliin 
the western part of the Twin City 7-County Metropolitan Area (Figure 1) . 
The entire Metropolitan Area ha i:)een mapx)ed previously through airphoto 
interpretation by tire Remote Sensing Laboratory under contract with a con- 
sortium of funding agencies.^ (In addition, in-house funds had already 
been used to perform tire single-date Landsat cinalysis of the study area.)^* 
The large eimount of existing data on this area provided a good test base 
for the study. The test site contains one of tlie most concentrated areas 
of wetland diversity in the v;est half of the Metropolitan 2\rca. It is a 

^The Himiesota Depcirtment of Natural Resources, the u.S, Army Corps of 
Engineers, the U.S. Fish and VUldlife Service, the USDA Soil Conserva- 
tion Service, and the ^Tv/in Cities Metropolitan Council. (VJerth et al . , 
1977; Owens and Meyer, 1978). 

^The University of Minnesota Agricultural lt<periment Station and the 
College of Forestry. 
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Figure 1. Study Site (loc^itod in the wootorn part of Hennopin County, Minnosota) 
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complex, heterogeneous area of undulating glacial topography reprosenta- 
time of a transition zone between tall grass prairie and northern hard- 
wood forest. The landsca]pe is interspersed with small farms, villages, 
cities, subdivisions, parks and golf courses. 

As the tests of tne various LANDSAT data analysis techniques were 
being apx^lied to the f> quadrangle test area, a para.llel test of digitised 
color infrared photographic data v;as undertaken in a limited svib-area (see 
Figure 1) . The digitized photographic data virtually removed the spatial 
resolution restrictions inherent in the LANDSAT data. These restrictions 
were found to be acute due to tlie spatial complexity of land cover in 
this area. In this way, the XDhotographic analysis permitted ci more basic 
test of the spectral characteristics and separability of the wetland 
types. This was deemed important for judging the potential applicability 
of Landsat data analysis procedures in areas vrhich are less complex spa- 
tially than those in tlie Metropolitan study area. It also provided some 
insight into the potential advantages to be realised with the improved 
spatial resolution of the Thematic Mapper on LANDSAT-D. In addition, the 
photographic research permitted cin initial investigation into the poten- 
tial utility of using digitised color infrared photography per se as a 
data source for v;etland mapping in spatially complex areas. 

The methods and results of the LAl^DSAT test and the test of the dig- 
itized photography comprise, respectively, the next two sections of this 
report. 
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IiANDSAT DATA ANALYSIS 

As previously indicated, loth single-date and double-date LANDSAT 
data sets v;ere analyzed in this study. Both data sets were processed 
at the Purdue University Laboratory for Applications of Remote Sensing 
(LARS). The single-date analysis was performed on data acquired on 
July 6, 1976. A coarse geometric correction (i.e., not using ground 
control points) v/as applied to this data set. The correction resulted 
in an average pixel positioning error of approximately 2 pixels N-S and 
4.5 pixels E-W v;hen image data v?ere compared to a 1:24,000 map base. 

This correction uses a transformation model based on nominal satellite 
parameters cind the geograiphic position of the test site, to rescale, 
deskew and rotate the original data. Following the geometric correc- 
tion, classification was performed using the IjARS hybrid multi-block 
cluster classification procedure (Hoffer and Fleraing, 1978) . This pro- 
cedure entails selecting blocks of training areas containing heterogen- 
eous cover types, clustering tliese blocks individually (unsupervised) , 
identifying the informational identity of cluster classes, pooling train- 
ing statistics frora similar classes in various blocks, and using these 
statistics to perform a maximum likelihood classification of the entire 
area of interest. Following this procedure, thirty-five spectral clas- 
ses were extracted from the single-date data to represent the following 
information classes: Non-v/etland, Forest, Shrub, Emergent, Submergent 

and Water. Overall classification accuracy as measured from 2555 ran- 
doraly selected x:>ixels of verified identity v/as 72%.. (That is, the total 
number of pixels classified correctly, 2555 = 0.72). However, tl-ie 
individual classes had highly variable classification accuracies. For 
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exaitiplC/ 96% of the Water test pixels were classified correctly, whereas 
only 10% of pixels known to be of the Shrub class wore classified cor- 
rectly. The average of all of the class accuracies v;as 41%, which v;as 
deemed inadequate for the application at hand. 

The generally poor classification perforinance and geometric correc- 
tion results of the single-date analysis v;ere improved U£oon in tlie double- 
date analysis process. First, application of the LARS geometric correc- 
tion, using ground control points, resulted in average additional pixel 
positional corrections of approximately 0.5 pixels M-S and 0.75 pixels 
E-W. The geometrically corrected data for twx) dates v/ere registered and 
the multi- block cluster tecliniqiie v?as used for classification. (In this 
case the July 5, 1976 data v/ere merged v/ith those of August 7, 1975) . 

The LARSYS feature selection process was employed to judge v/hich of the 
eight available bands of data was optimum for subsequent classification. 
Using this process, bands 5 and 6 from both dates were deemed best suited 
for the classification. Given the jpoor performance of the Slirub class in 
the single-date analysis, this class v/as eliminated from the double-date 
classification. Also, it was hypothesized that separation between the 
Hon-wetland class (basically corn) and Emergents (cattail) v/ould improve 
with the introduction of the August data, since the corn had tasseled by 
this date. 

The test field accurcicy of the double-date classification v/as eval- 
uated in a manner similar to the single-date effort. This involved eval- 
uation of some 450 randomly selected test fields that varied in size from 
one to nine homogeneous pixels. The cover type present in each field v/as 
verified by ground site visits. Table 1 summarizes the results of com- 
paring the LAMDSAT and ground data. As can be seen from Table 1, the 
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Table 1. Test: Field Accuracy Assessment of Double-Date LANDSAT 
Glassificationt 


Knovm Cover lypes 

% Clas 

sified As 

Indicated 

Cover Type 

Cover Type 

No. of 
Test Pixels 

cS 

H 

■P 

1 

Forest 

Emergent 

u 

QJ 

4-» 

Nonwetland 

1,706 

75 

8 

17 


Forest 

120 

34 

56 

10 


Emergent 

121 

33 

30 

35 

2 

V7ater 

523 



2 

98 


Overall Classification Accuracy: 77% 

Average Accuracy Per Class: 66% 


tA detailed discussion of the classification and mapping accuracy of the 
single and double-date LANDSAT cinalyses is given in V7erth (1930) . 
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overall classification accuracy for the double-date analysis was 77 ^j. 

The accuracy for the various information classes ranged from a low of 
35% for tl'ie Emergent class to 98% for Water. The per class average accu- 
racy was 66%. The comparative levels of errors of omission and commis- 
sion for each class are also shov/n in Talkie 1. 

In short, the two-dcite results indicated an increase in overall 
classification accuracy of approximately 5% (from 72 to 77%) and an 
improvement in per class classification accuracy of about 25% (from 41 
to 66%) . Additional classification accuracy appears to be precluded by 
the spatial and spectral resolution limits of the LAMDSAT MSS relative 
to the complexity of the study eirea. This is not to say that LANDSAT 
has limited application to v;etland mapping in general. Rather, we have 
simply concluded that one should expect poor success in single-date anal- 
yses and only moderate success in double-date analyses of vretland areas 
v/hich are as complex as those in our study area. 

ANALYSIS OF DIGITIZED PHOTOGRAPHY 


As mentioned previously, digitized color infrared aerial photography 
was also analyzed in tliis study as an extension of the LANDSAT analysis 
effort. Again, this source of digital data virtually removes the spatial 
resolution restrictions limiting tlie LAITDSAT automated classification 
process, thereby permitting more tliorough study of the spectral proper- 
ties of the wetland cover types under investigation. Constraints of time 
and funding limited the photographic analysis to coverage of a single 
study site located \vithiri the LANDSAT test area (Figure 1) . This site 


/Y 


f 


AID 


contains ci nuinber of small v;etland areas v/hich generally characterize 
those of the surrounding region. 

Figure 2a is a (1.7>:) enlargement of the original 70mm £>osition 
transparency used in tlin photographic cinalysis. The original image v;as 
a 1:50,000 photograph (Kodak type 2443) taken writh a Hasselblad 500 EL 
(f=50mm) CcUi\era on July 21, 1977, in support of previous wetland research 
efforts conducted by the Remote Sensing Laboratory. The image was digi- 
tized using a modified P-1700 Optronics drum scanner made available by 
the University of Wisconsin Environmental Monitoring and Data Acquisition 
Group (UV7/EMDAG) . Density readings from 0-3D v;ere digitized into 256 
levels and recorded on computer compatible tape. Broad band separation 
filters were used to obtain the approximate image density record for each 
of the three film layers. A measurement spot size of 50ym v;as employed, 
so that each pixel represented a ground area approximately 2.5m square. 
The spectral and spatial quality of the resulting digital data can be 
seen in the color rendering of the digitized data shown in Figure 2b. 

This figure was generated with the University of Minnesota D-47 Dicomed 
Color Image Recorder. The digital data from each image band were con- 
trast enhanced by a histogram equilization algorithm prior to being 
recorded. This algorithm alters the distribution of original image 
values so that each display level comprises approximately an equal area 
portion of the disx:>lay iinage in each band of data. Figure 2b represents 
a digitally generated color composite of the three original bands 
enhanced in this manner. 

The photographic analysis effort took on tv;o forms: 1) a full- 

scene classification, the intent of which was to classify all pixels 


and cover types present in the scene, and 2) a sub- scene c l assification 



to clar.Gify v;ctland area?; and cover types only. In both catios, a super- 
vicucl classification approach was cnploycd using the *Tppior Oi^timiza- 
tion*’ fc>"m of the naximum likelihood classifier that is inte.gratod into 
the University of Minnesota Image Processing Software (UfllPS) . A dos- 
crixjtion of the two classification procedures follows. 

Full-Scone Classification 

Training for the photographic classifications was accoipplished by 
photo-interpreting tx*aining areas on prints produced from the cnlianced 
Dicomod image. Tlie boiuidary of each training area \;as then digitized in 
terms of x and y image coordinates. In turn, these coordinates w’ere con- 
verted to coliurn and row addresses in tlie scanning microdonsitome tor data 
using a i:>olynonial two-dimensional coordinate transformation. This 
transformation took the form 

columji no.«a +ax*fay + a>:^ + aY^ 

12 3 5 6 

row no. « a + a x -f a y + a xy + a x^ + a y^ 

7 8 9 10 11 12 

A least squares observation equation solution for coefficients 

a ...a was performed by measuring both the image and line printer out- 
2 12 - 

put coordinates of readily identifiable ground control x'^oints in the 
scene. Using 10 xx^ints for this x^urj^ose, tlie root meeui square error 
(rms) of the transfonration v;as found to be 0.89 pixels in the x and 0.74 
pixels in the y image direction. 

Training statistics v.ero dovclcx^cd for the 13 miajor cover typos 
ax>X'V :ring in the scene: Uator, Duckweed ( Lemna minor ) , Water Lily (Uym- 
phaoa spx->.- IJiix^har spp. ) , Cattail (T>pha spp.), Reed Canary Grass ( Phal- 
aris arundinacca ) , Vlillow ( Salix spp. ) , Forest, Cropland, Pasture, Bare 


Soil, Road, Roof, and r*hailow. Multiple npoctral clacuos (25) were de- 
rived froir. 35 trainiriy iirodB to ad:^quatciy train for the 13 cover typoD. 

Tabic 2 iu a couf union ratrix indicatinq the training field clacs- 
ification accuracy obtained in the full-nceno analysis. Note that for 
BOLie tyx^es the accuracy of classif icaion in tlic training areas v/as excel- 
lent (e.g., Pasture « 97^w) . Other types, most notably \7illow, were very 
poorly classified (25%) . Figure 3 is a color-coded Dicomed image sliow- 
ing the results of apx>lying the full-scene classification. As can be 
seen from both Tabic 2 and Figure 3, similar to the case with the LANOr.AT 
data, many of tlic cover tyx^es i^rcuont in tlie scene could not be sx>ect- 
rally sex>aratcd. However, of particular note is the fact tliat virtually 
all wetland types wore mutually separable. Tliat is, errors of ondssion 
and coinnission in each wetland type v;ere primarily attributable to the 
spectral confusion between v;etland and ncn-v/etland cover tyx'os, rather 
than betwec»n wetland cover types themselves. 

Sub-Scene Classification 

Because of t)ie spectral confusion betv/een the wetland and non-wet- 
land cover types, it v;as hyp>othesizod that if the training £ind classifi- 
cation procedures v7ore limited to v/etland areas only, the classification 
accuracy w’ould improve greatly. Hence, the training process was repeated 
in wetland cireas only, using tlie following classes: Duckweed, V/ator, 

\/ater Lily, Cattail, Heed C^mary Crass, and v;illow. The boundaries of 
the wetland eireas w^erc delineated again using conventional image inter- 
pretation. This process recjuircs much lest: labor than discriminating all 
cover tyx^e Loui:darios occurring within the wetlands. The pur^jose here \;as 
siiiply to separate wetland and i:on-wctland image areas. Sx^ectral pattern 


Vrbie 2. Claysif ication PorforiiViiicc for Training i’iolrls Usod in rull-i:cone 
Analyni s 


y'jiovm Cover Tyx>e 


t Clascified At; Indicated Cover *I'i'i)e 


>1 

u 


Cover Type 

No. of 

Train. Pixels 

I 

1 Duckweed 

1 

VI 

o 

H 

•H 

^ rH 

•H 

Vi rj 

U 4^ 

4J 44 

/J 

U U 

rJ 

c 

rJ 

U 

'O 

Cl 

o 

o 

.H 

rH 

;H 

j Forest 

1 

j Cropland 

Pasture 
Pare Soil 

^0 

r2 

Duckweed 

C7 

87 


2 




4 

1 

G 

\:ator 

1,700 


81 






1 


VJater Lily 

240 



89 




10 

1 


Cattail 

319 



91 


9 





Peed Canary 

C64 




74 


10 

If. 



Willow 

255 



17 

3 

25 


55 



Fr»rest 

48 



2 

4 


73 

17 

4 


Cl ^^vland 

3,392 



4 

2 

6 

1 

86 

1 


Pasture 

CIS 



1 




, 2 

97 


bare Soil 

232 

7 






11 

81 

1 

Road 

400 

2 






1 

3 1 

93 

Hoof 

5G 










Sliadow 

86 


7 









o 


18 


100 


93 


Overall Classification Accuracy: 
Average Accuracy Per Claes: 82 t 


84 ?. 


Shadow 



Color Ko 


Duckweed - magenta 


Forest - black 


Water - blue 


Cropland - tan 


Water Lily - magenta 


Pasture - green 


Bare Soil - brown 


Road/Koof - white 


Reed Canary - red 


Shadow - black 


Figu>e 3. Full-scone classif icati on of digital photographic data 



Table 3. Clacaif icatioa Performance lor Training Ficlcli; Used in 


Known 

Cover IVP^ 


Classified 

As 

Indicated 

Cover 

Type 

Cover Type 

No. of 

Train. Pixels 

X 

E 

Vi 

Q) 

vratcr Lily 

Cattail 

>1 

i 

•g 

V.'illow 

Duck\;e(^ 

67 

00 


4 

0 



Water 

1,700 


100 





V/ater Lily 

240 


. 

100 



Cattail 

119 

1 



89 


10 

Reed Canary 

664 





97 

3 

Willow 

255 




40 

21 

39 


Overall Claaaification Accuracy: 
Average Accuracy Per Clasa: 05.5^4 
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TablG 4. ciasEjificatiqn Porforitiauca in Eight VAitlancl Teat Areas 


Kno\'m 

Cover Type 

Z Clas 

sified A 

s Indicated 

cover Typo 

Cover Type 

No. of 
Tost Pixels 

rd 

0 

rQ 

u 

E 

M 

0 

$ 

u 

.M 

1 Cattail 

Eeed Canary 

Duclwced 

984 

95 


3 

2 


Water 

6,173 

2 

97 



1 

V/ater Lily 

2,240 

11 


72 

1 

16 

Cattail 

13,586 




81 

19 

Reed Canary 

9,920 

1 


5 

2 

92 

Overall ClasE 

lification A.ccuracy; 

87 ^ 






Average Accuracy Per Class: 87*!! 


recognition v/aa then used to provide the more aubtlo delineation of. wet- 
Icind species, a more time consuming process for the i;ho to- interpreter. 
Table 3 indicates the classification fccuraoy obtained in the "wetland 
only" training areas. Note the improvement in training field classifi- 
cation accuracy v/hi ;h results v/hen non-wetland classes are eliminated 
from the ar/alysis. Figure 4 shows the results of the sub-scene classi- 
fication. ,i,,, As can be seen from Tdsle 3 and Figure 4 the Willow typo 
introduced the largest errors in the suL-scene classification. 

T i a final classification effort, the Willov/ type v;as dropijed and 
the resulting statistics were again applied to the "wetland only" subset 
(with predominemtly v/illow areas excluded) of the original imago. The 
accuracy of this classification was evaluated in eight tost areas 
located in the scone. The result of the classification in these areas 
is shovm in Figure 5. The classification accuracy in each of the eight 
test areas v/as obtained by comparing the computer-derived classification 
to that obtained by photointerpretation. This was done for all pixels in 
the eight areas (approximately 33,000). In those cases v/here the photo- 
interpretation v?as ambiguous, a field visit to the site v/as made to be 
certain of the correct identification of the cover type present. 

The results of the pixel-by-pixel accuracy assessment are given in 
Table 4. Errors in tlie classification of Ductaveed and Water are be- 
lieved to be caused by edge effects, in v/hich a given pixel actually 
overlays tv/o cover types near a boundary between the two types. To 
some extent, this effect enters iiito the accuracy problems for Water 
Lily as v/ell. It 4s felt that the classification accuracy of the Water 
Lily and Reed Canary classes could be improved som.ewhat with additional 
training. Most of the Water Lily/Reed Canary confusion occurred in an 


Color Key 


iquro 4 


Duckweed - magenta 
VJatcr - blue 
Water Lily - magenta 
Cattail - orange 
R<'‘ed Canary - red 
Willow - yellow 

Sub -scene classification of photgraphic data in v;etland arc^a; 
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area of tlic nceno which wan not used for trainirej. The errors in class- 
ification of Cattail arc believed to be due to the presence of nixed 
Cattail/h’illow in irany areas. A sparse scattering of V/illow occurred 
in rany areas of both the Cattail ai.c tlic Rood Cai\ary cover types. 
\7licrcvc:r present, V/illov; tended to cause nisclassification due to the 
highJy v^ariable spectral characteristics* of this cover tyx>e. These 
difficulties may be resolvablci the inclusion of a ’’texture” variable 
in future anaJ.ysis i^rocedures. 

Based on the results of tliis investigation, the follovv’ing general 
conciusions have been reached: 

1. If, within the accuracy levels it is capable of j^rovidiiig, L?.NDdAT 
data are judged to be accci>l:able for v/etlands classification, a 
precision geometric correction should be ax)plied and multi-date 
data sets should be used. In this study a 25% improvement in aver- 
age classification accuracy was realized by x'^^^ocossing double-date 
vs. single-date data (from 41% to 6G%) . Under the spectrally and 
spatially complex site conditions characterizing the geographical 
area used in tliis study, further iiaprovem.ent in w'otland classifi- 
cation accuracy is ap;>arcntly precluded by the sx^octral raid spatial 
resolution restrictions of the LANDSAT MFf». 

2. Full-scene cUialysis of scanning donsitometor data extracted from 
snail scale color infrared xjhctograx^jhy failed to i^emit discrimi- 
natioii of many w^i'tland and non-wetland cover tyx:?es. This v/as due 
cssentjally to tlie .Seine sx^ectral confusion yxitwecu cover tyx-es that 
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was realised in the LANDSAT analysis. When classification of photo- 
graphic dab. was limited to wetland areas only, much more detailed 
and ciccurate classification could be iricide. The final subset class- 
ification made in tliis study yielded an average classification 
accuracy of 87% witih a 2.5m spatial resolution. 

3. The integrcition of conventional iraage interpretation (to simply 

delineate v/etland boundaries) and machine-assisted classification 
(to discriminate among cover types present within the wetland areas) 
appears to warrant further research. Though both the delineation 
of wetland bovmdaries and tlie development of valid training area 
statistics are labor-intensive tasks, the results from the combined 
approach appear to be much more accurate and economical than those 
obtained in conventional full- scene digital cinalysis. Additional 
research is needed to study the feasibility and cost of extending 
this methodology over a Icirge circa using LANDSAT and/or small scale 
photography. 
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